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Abstract. In this paper we discuss a new strategy for the detection of gravitational 
radiation likely emitted by cosmological gamma ray burst. Robust and conservative 
estimates lead to the conclusion that the uncorrelated superimposition of bursts 
of gravitational waves can be detected by interferometric detectors like VIRGO or 
LIGO. The expected signal is predicted to carry two very distinctive signatures: the 
cosmological dipole anisotropy and a characteristic time scale in the auto correlation 
spectrum, which might be exploited, perhaps with ad hoc modifications and/or 
upgrading of the planned experiments, to confirm the non-instrumental origin of the 
signal. 

PACS numbers: 04.80.Nn, 95.55.Ym, 98.70.Rz, 98.70.Vc 
1. Introduction 

Gamma Ray Bursts (GRB) are short and intense flashes of e.m. radiation in the 
keV-MeV range, accidentally discovered more then 30 years ago 0. The observation 
with high spatial resolution by the Italian-Dutch satellite Beppo-SAX || of the X- 
ray afterglow, has allowed the optical identification of some of the burst, which has 
given the evidence of the cosmological origin of the GRB's. It is now evident that 
for few seconds, at least twice a day, unknown sources become so bright to over 
shine all the visible Universe [[J. From the observed 7-ray fluence and the red shift 
of the host galaxy, it is estimated H that the average buster emits isotropically 
{E liso ) = (l.3tj' 2 ) x 10 53 erg (We will assume |J here and in the rest of the paper 
a flat Universe with Hq = 65 km s _1 Mpc -1 , Qm = 0.3 and Qy = 0.7). The dispersion 
around this average is consistent with a log-normal distribution having logarithmic 
width a~ jiso = l-7llo 3 (about a factor 50). This energy corresponds in the average to 
« O.lM c 2 , but arrives in the extreme case of GRB 990123 to be ~ 1.4 M c 2 , which 
poses serious problems for understanding what could be the "central engine" powering 
these events. In the currently accepted "fireball" model JI| the 7-rays are emitted by 
electron/positron pairs, accelerated by internal relativistic shocks, which radiates via 
synchrotron and/or synchro-Compton process. In order to be optically thin, the fireball 
should expand with ultra-relativistic velocity (bulk Lorentz factor T »s 100). Even if 
the conversion of the kinetic energy into 7 radiation is estimated f| to be very efficient 
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[rj 1 pa 0.2), the formation of the relativistic wind with a kinetic energy pa Eyiso/rj^ could 
be problematic. 

The energy output of the central engine of GRB can be largely reduced if the 7-ray 
emission is beamed with a small opening angle. This possibility is supported by the 
observation of an achromatic break in the time evolution of the afterglow luminosity, 
occurring when the bulk Lorentz factor becomes T < 9j^ t . In a recent paper || this 
hypothesis is applied to estimate the opening angle of the 10 GRB's with observed 
afterglows, with the interesting result that the average opening angle in the sample is 
(9 j e t) — 4°. Furthermore, it appears from these data that the fluence of the GRB in the 
sample is correlated with the opening angle, in a way that justifies the large spread of 
the apparent fluence as entirely due to the spread of the opening angles. The intrinsic 
energy emitted by GRB in 7-rays, if estimated as 

E ljet = f jGt F 7 d 2 L (l + z)- 1 (1) 

where fj et = \ (1 — cos^t) is the beaming factor 0, clusters around a value of 
(Eyj et ) = 5 x 10 50 erg with a logarithmic width o" 7jet ~ 0.3 that is about a factor 
25 smaller than the one of the E^i SO , If this is true the kinetic energy would be of the 
order of Eqrb ~ 1O~ 3 M c 2 , about constant for all the bursts. The scenario emerging 
from this study is particularly appealing because one could assume that the central 
engine of the GRB produces always a similar amount of energy, while the complexity of 
the energy transfer from the central engine to the relativistic jet causes the wide range 
of opening angle, and consequently the wide range of apparent fluence that is observed. 

This scenario implies also that the true rate of bursts Rj et ~ (fjet) Riso, is enhanced 
by a factor (fj e \) ~ 500 respect to the observed rate. However clearly the total 
energy output per unit volume and unit time injected in the universe by GRB explosion 
e 7 = Rjet E 1 jet = RisoE^i SO will be the same of the isotropic case. 

In the next Sect. |^ of this paper we discuss the implications of the possible beamed 
emission on the detectability of gravitational waves (GW) likely emitted by the collapse 
event that originates the GRB. It is to be expected that if the e.m. energy of the 7- 
ray burst is not extraordinary (about the same order of magnitude of the total output 
of a Type lb SN), the energy output in gravitational waves should also be modest. 
But, as we observed above, in the latter case the number of bursts should be increased 
correspondingly. This suggests a new strategy for looking to gravitational radiation from 
GRB. As we will show in Sect. |3|, even if the GW's produced by each GRB are below 
the detection threshold, integrating over one year one should find an excess of noise, 
due to the uncorrelated superimposition of many GW pulse trains. In Sect. |]we discuss 
two very distinctive signature that a genuine cosmological stochastic signal should carry. 
Those intrinsic signatures, perhaps detectable with modifications and/or upgrading of 
the planned experiments, could be of great help to check the non-instrumental origin of 
the stochastic signal. Finally in Sect. |5| we summarize the results of this investigation. 
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2. Gravitational wave signal from GRB 

The energy required to power the 7-ray burst, estimated by Eq. ([[]), is of the order of 
« 10~ 3 Mq c 2 , not overwhelming but always much larger then the total nuclear binding 
energy of few-M 's stars. Therefore it requires in any case an energy release that is 
compatible only with the collapse of a several solar masses object. A natural assumption 
would be that the GRB are similar to ordinary supernovas, but perhaps the difference is 
in the final result of the collapse. Like for the "failed supernova" model fl6| the collapse 



of a massive star to a black hole surrounded by a dense rotating torus of material that 
might result in a relativistic jet. The energy irradiated in this case as gravitational 
waves can be estimated [jDJ as E GW rs e GW M BH c 2 where M BH is the mass of the black 
hole and e G w ^ 10~ 4 is an efficiency parameter. As we observed above the bolometric 
energy released in the gamma ray burst is already of the order of E^j et > 10~ 4 M & c 2 
and 7] 1 0.2, then we could expect that the energy provided by the central engine 
could be Egrb ~ £grb Mbh c 2 where €grb is the fraction of the gravitational energy 
converted into kinetic energy of the ejecta subsequently radiated as 7-rays. 

It could be possible under certain conditions |18[ that €grb ~ £gw, but perhaps 
we should adopt the more conservative view that e G w ^ £grb even if it is unlikely that 
e GW <S ^grb- Our ignorance of the actual mechanism that could transfer the energy 
from the accretion torus to the jet does not allow a solid prediction of the energetic of 
the central engine of the GRB. Nevertheless we can set a lower limit to this energy in 
the form Eqrb <i E 7 j et /r) 7 . We can parameterize the production of GW introducing a 
phenomenological r/ GW < 1 putting E GW = rj GW E^j et /r]^ and verifying the sensitivity 
of the GW detectors to the actual value of r\ G w- 

The energy flux (viz. energy per unit surface and unit time) of GW produced by 
a burst of intrinsic luminosity (in gravitational waves) L G y/ at a red shift z is given by 

El 

b (1 + z) L GW 

Fgw ~ And 2 L (2) 
where d^ is the luminosity distance. Assuming a typical time scale for the emission At 
we can estimate this flux at Earth for a typical red shift z = 1 

fjet E^i s0 / Tj^ 

2.5 x 10 51 eig / 



F GW {z = 1) » 4 x 10- 7 f j'f^I ' n ) f^) (^] ergcm-VSr- 1 



fVGw\ ( 10 ms\ 

— : erg < 

V0.01/ \ At J 6 

In order to convert the flux of Eq. (0) into an adimensional amplitude we can use the 



classical formula 20 



c 3 



hl + hl) (3) 



16nG 

where the average is taken over several wavelengths. The amplitude of the signal 
produced depends from the direction and the beam pattern of the detector. In the 
best case we have integrating over time and applying the Parseval' s theorem 

" 2~ h 2 ( \ j 16ttG (1 + z) E GW 
uh b (u J )du = — (4) 
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where hl(uj) will be the Rayleigh power of the signal as a function of the frequency uj. In 
order to estimate the order of magnitude of the amplitude of the GW signal we do not 
need a detailed shape of the spectral power density of the signal, but only the knowledge 
of the firsts two moments of the distribution uj and Auj. In fact we can recast Eq. (|4]) 
in the form 

16vrG (1 + z) E GW 



/+oo 
h 2 b (u)du 
-oo 



C 3 47T d\ 



It is remarkable that rather natural physical assumption on the first and second moment 
of the unknown Rayleigh power distribution h\[uS) can be made. In fact we can assume 
that the first moment will be uj ~ 2nc/rs where rs will be the Schwarzschild radius of 
the collapsed object The second moment will be the r.m.s. bandwidth of a wave packet 
that can be estimated Auj « 2k /At where At is the duration of the emission in the 
comoving frame. In this case we have, if cAt 3> rs, for the peak amplitude 



The expected amplitude for typical values is 



rf^r-lO- 27 (^-) ( EgR Z V (-**-)* Gpc/VKz 

b V5 km/ V2-5 x 10 51 erg J VO.Ol/ V 10ms / ' 

It is clear from this estimate that the probability of detecting a single burst is very low, 
unless t]gw 3> 1- However even if the individual burst could not be detected, it is to be 
remarked that in case of beaming the rate of explosion would be very large (~ 1500 per 
day), therefore the stochastic accumulation of signal integrated over a long time could 
emerge from the noise. 

3. Cosmological background 

The energy flux of GW produced at Earth from a cosmological distribution of sources 
is given by 

**' = jf R {z) E OWJT - m7) (6) 

where R(z) is the comoving rate of bursts per unit volume exploding at red shift z, E GW 
the average energy emitted in gravitational waves by each source that we have estimated 
in §0 and H{z) = R/R is the Hubble expansion rate fT9[| . As obvious, the lower apparent 
brightness of distant burst is compensated by the increase of volume up to a red shift 
1-2, like in the classical Olbert' s paradox. It is also worth noticing that the flux will 
be the same for beamed or isotropic sources. The rate R{z) can be obtained from the 



Log N-Log P distribution |T(J, |llfl , even if its value depends from the assumption made 
on the evolution of the burst rate in the recent past. In practice the energy flux of the 
present stochastic background of GW is dominated by the rate of explosions at z > 1. If 
the GRB are produced by a final catastrophic collapse event of an evolved massive star, 
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one expects that the large scale distribution of GRB should reflect the star formation 



rate. The latter is known |12], [13| to increase by at least a factor ten from z = to 
z ~ 2 to flatten in the range 2 < z < 3 and to decrease exponentially for z > 3, reaching 
a value similar to the present for z ~ 6. The fit obtained with this evolution is 
Riso(z = 0) = 0.14 ± 0.02 Gpc _3 y _1 . On the other side if no evolution is assumed the 
constant rate is about one order of magnitude larger. But it is worth noticing that due 
to the constraint of having about 1000 bursts per year the only difference between the 
two cases is that in the first case the average distance of the bursts is slightly larger, due 
to the fast rise of the rate of explosion with red shift. In fact from Eq. (P) we estimate 
for constant rate 

F dlfS « 10- 10 ( Eiso/f]l ^ (^\ ( — ^ ere cm" 2 s- 1 sr" 1 

* GW iU UsxlO^ergJ WHlGpc-VV 

while assuming evolution of the rate in the recent past (see below) the predicted flux 
is practically the same. Comparing with Eq. (0) we observe that the diffuse flux is 
consistent with the order of magnitude estimate 

F%w/ F gw « (fiel) x 1,000 bursts/1 year x At 



Applying again Eq. (0) we have: 



-.3 



dz 



h\ + h\ ) = / R(z) E GW (7) 

On the L.H.S. of this equation we have the amplitude of the wave that invests at a 
certain instant the detector. We have seen in the previous section that each of this 
burst will not have a detectable intensity, but if we average over an observation time 
T long compared to the GW burst duration but short compared to Hubble time scale 
(typically one year) we have a signal 

that will be detectable if the power spectral density is greater then the power spectral of 
the noise, averaged over the same observation time. The uncorrelated superimposition of 
bursts of gravitational waves will be well approximated, for the central limit theorem, by 
the superimposition of red shifted gaussian distributions. Therefore the power spectral 
density of the signal can be estimated by the integral 

,72, x\ G f ZmaX r>, ^ E GWr 2 s l(i+-)"-"ol a dz 

(h^uj^T^^— / R(z) t= — e ^ — — (9) 

^ K " vr 2 c 5 in 1 ^v^Aa; (l + z)H(z) KJ 



o 

where the normalization has been obtained from Eq. (^). We have reported in Figure 
which is the central result of this paper, the Rayleigh power of the stochastic signal from 
the cosmological background for tjgw = 0.01. The dashed curve in Figure |I] represents 
the power expected allowing for an evolution of the rate R(z) with the red shift similar 
to the luminosity density evolution of QSO's This evolution is assumed to be 



R(z) = (1 + z) a with a ~ 3 at low red shift, z < 1.9, a = for 1.9 < z < 2.7 and 
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Figure 1. Stochastic signal integrated over one year from cosmological GRB 
(vgw = 0.01). In this plot the solid line corresponds to a flat rate, while the dashed 
one to a rate of bursts evolving like the QSO rate (see text). The dot-dashed line is 
the noise estimated in the VIRGO proposal pa], averaged over one year. 



an exponential decay at z > 2.7 [Q, similar to that describing the evolution of star 
formation rate that we have discussed above. The solid curve represents, on the contrary, 
the power expected if the rate of the GRB is assumed to be constant. For comparison 
we have also reported the noise expected in the VIRGO experiment [25] averaged over 
one year of integration. This comparison shows that the cosmological signal should 
be detectable at the same level in both cases, because the overall normalization is 
constrained by the observed rate by BATSE of ~ 1000 bursts per year. 



4. Signatures of the genuine signal 

We expect for the stochastic signal produced by cosmological sources some very clear 
signatures showing its origin. We will not discuss in the following whether or not 
those signatures will be detectable by the two planned detectors VIRGO and LIGO, 
but we expect that an ad hoc modification of those experiments could be designed in 
order to extract all the possible informations from the signal. The more evident of 
those signatures should be the dipole anisotropy, as observed for the cosmic microwave 
background (CMB). The dipole anisotropy of the CMB is interpreted as the result of 
Doppler shift caused by the solar system motion relative to the isotropic radiation field. 
This motion is confirmed by measurements of the apparent velocity of local galaxies 
p6f . The motion of the observer (receiver) with velocity (3 = v/c relative to a source 




Figure 2. Amplitude of the dipole anisotropy intrinsic modulation on the GW 
background as a function of the frequency. Solid line is for GRB rate evolution, while 
dashed line is for constant rate. In this plot a negative amplitude corresponds to a 
modulation in opposition of phase. 



of gravitational waves with frequency uq produces a shift in the observed frequency u' 



given by the formula |27 



= , (1 P\ « uj [1 - (3 cos 9 + 0((3 2 )] (10) 

1 — p COS U 

where the velocity for the solar-system barycentre is |28| , |29[| (3 = 0.001237 ± 0.000002 
or v — 371 ± 0.5 kins -1 and 9 is the angle formed with the direction of equatorial 
coordinates {a, 5) = (11.20 h ± 0.01 h , -7.22° ± 0.08°). This frequency shift produces at 
a given frequency a diurnal modulation of the Rayleigh power in the true signal only, 
which depends from the slope of the spectrum. From the Figure p] we can infer that this 
modulation will be particularly enhanced at the extremes of the range of cosmological 
red shifted distribution of the characteristic spectrum. In Figure [| we have reported 
the theoretical maximum amplitude of the intrinsic modulation, induced by the dipole 
anisotropy as a function of the frequency, calculated substituting u' Q given by Eq. ([10|) 
to u into Eq. (H). 

A second signature of the cosmological signal comes from the fact that the 
detectable power is obtained accumulating many short duration GW pulse trains. The 
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signal auto correlation spectrum denned as 

A(t) = i J\{t) h(t + t) dt (ll) 

should show the evidence of a correlation over the characteristic scale of the stellar 
collapse (of the order of 10 ms.) If the instrumental noise is white and the sampling 
frequency of the detector very high, the distinction could be very clear. In real life the 
detection of this feature could be much harder because the noise of the detector will 
be the superimposition of white noise due to microscopic processes and colored noise 
coming from other physical sources (like for example the mirror resonance or the 1/f 
noise) and the sampling frequency could be inadequate. 

5. Conclusions 

In §Q we have shown that the GW emission from single bursts at cosmological distances, 
if the e.m. emission is beamed with a small angle as suggested by afterglow observations, 
is expected to be, for conservative values of the source emissivity, d L h < 10~ 27 Gpc/v^Hz 
which is well below the detection threshold of presently planned experiments. However, 
as we have discussed in §|3|, if the e.m. emission is beamed only a small fraction (one over 
500) of the GRB are observed by 7-ray satellites. This implies that small amplitude 
pulse trains of GW impinge over the detector at a frequency that could be as high as 
one per minute, which is practically a continuous signal. Even if the individual pulse 
is small, integrating over a reasonable observation time (order of one year) an excess 
Rayleigh power should emerge from the instrumental noise. The frequency spectrum 
of this eventual excess should give a direct information on the characteristic time scale 
of the collapse and on the cosmological evolution of the GRB rate in the recent past 
(z « 1 — 2). The predicted amplitude is conservatively estimated to be in the range of 
5 x 10~ 26 1/vHz averaging the Rayleigh power over one year. This estimate is rather 
robust because does not depends on the beaming factor and depends only slightly from 
the large scale distribution of the GRB sources. The cosmological origin of the excess 
noise can be proved by detecting a dipole anisotropy, that in the relevant frequency band 
could reach the level of a fraction of percent. In addition the auto correlation spectrum 
of the noise shall carry an imprint of the characteristic duration of the GW pulse trains. 
The detection of those two signatures, even if perhaps not possible with presently 
planned experiments, can give the important additional evidence of the cosmological 
origin of the stochastic signal and informations on the physics of GRB's. 
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